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We discuss, from a theoretical point of view, the four 
wave mixing spectroscopy on an ensemble of p-doped 
quantum dots in a magnetic field slightly tilted from the 
in-plane configuration. We describe the system evolution 
in the density matrix formalism. In the limit of coherent 
ultrafast optical driving, we obtain analytical formulas 
for the single system dynamics and for the response of 
an inhomogeneously broadened ensemble. 



The results are compared to the previously studied time- 
resolved Kerr rotation spectroscopy on the same system. 
We show that the Kerr rotation and four wave mixing 
spectra yield complementary information on the spin dy- 
namics (precession and damping). 
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1 Introduction. The properties of confined spins in 
semiconductor nanostructures are interesting because of 
their expected important role in quantum computing and 
spintronic devices. In particular, the hole spin attracts much 
attention because of its enhanced coherence time. In sys- 
tems that efficiently couple to optical fields, like quantum 
wells and self-assembled quantum dots, an interesting pos- 
sibility is to study the spin dynamics using optical spec- 
troscopy tools. A method which is particularly suited for 
this purpose is the time-resolved Kerr (or Faraday) rota- 
tion (TRKR, TRFR) QHD, where the evolution of the 
occupations of the hole and trion Zeeman states is traced 
by investigating the rotation of the polarization plane of the 
reflected or transmitted probe pulse. 

In a recent work J4], we proposed a general descrip- 
tion for the dynamics of a confined hole-trion system in a 
tilted magnetic field. We showed that a TRKR experiment 
provides rich information about the rates of spin precession 
and decoherence (both longitudinal and transverse). In par- 
ticular, the optical response at slightly tilted magnetic fields 
contains contributions related both to the hole spin relax- 
ation (longitudinal decoherence) with respect to the hole 



spin quantization axis, as well as the dephasing of spin co- 
herences (transverse decoherence) with respect to this axis. 
Under favorable experimental conditions (a sufficient sep- 
aration of time scales), the corresponding two decoherence 
rates can be deduced from a single run of an experiment, 
although in a realistic system the latter may be convoluted 
with (or even dominated by) the inhomogeneous dephas- 
ing. 

Although the TRKR or TRFR method may be the most 
obvious choice for the investigation of the spin dynamics 
the spin precession of the trion and hole in a tilted field 
will lead to transitions between optically active and in- 
active states which should be manifested in any form of 
the optical response. In particular, signatures of the spin- 
related dynamics should be visible in the four wave mixing 
(FWM) nonlinear spectroscopy. As the FWM spectroscopy 
is one of the most widely used methods in the investigation 
of semiconductors and their nanostructures [5 6 1 it may be 
interesting to study how the spin precession and decoher- 
ence affect the FWM response and what information on 
the spin-related kinetics can be extracted from this kind of 
experiments. 
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In this contribution, we study the FWM response of the 
same system of confined holes as discussed in our previous 
work J4j. We show that the FWM signal is also affected 
by the spin dynamics in a way that allows one to extract 
the Larmor frequencies of electrons and holes. The avail- 
able information on the decoherence rates is less specific 
than that contained in the TRKR response and is available 
only as long as the hole spin dephasing times are not much 
longer than the life time of the optical coherence. On the 
other hand, the spins undergo partial refocusing in the two- 
pulse optical echo experiment which leads to the appear- 
ance of components in the optical response that depend on 
the inherent spin dephasing rates but are insensitive to the 
inhomogeneous distribution of the Larmor frequencies. In 
this way, the FWM method may be a valuable experimen- 
tal tool to study the spin dynamics, complementary to the 
TRKR or TRFR techniques. 

The paper is organized as follows. In Sec. [2] we define 
the model under study. In Sec. [3] we derive the FWM re- 
sponse of the system. Next, in Sec. [4] we discuss the result 
and compare the spin-related information contained in the 
FWM signal to that available from the TRKR response. 

2 Model. We study a system similar to that discussed 
in our previous work Q, consisting of an ensemble of 
quantum dots or trapping centers in quantum wells with 
one confined hole in each of the dots. The system is excited 
with a pump pulse at the time t = —t and a probe pulse at 
t = 0. Unlike in our previous work, we are now interested 
in the FWM response from the system. Experimentally, the 
relevant third-order response is isolated by choosing the 
appropriate excitation and detection directions (see Ref. J4) 
for a discussion). In the modeling, we assume the pump 
and probe pulses to have phases <f>\ and 02, respectively, 
and calculate the terms in the response of an inhomoge- 
neously broadened system that carry the phase 2^2 — <t>i ■ 
Out of many possible configurations of the polarizations of 
the excitation and detection, we choose the circularly co- 
polarized one, with both pulses having the cr+ polarizations 
and the detection being performed at the same polarization. 
The system is placed in a magnetic field tilted by the angle 
d to the normal, which defines the trion spin quantization. 
The hole spin is quantized along the axis at an angle tp from 
the normal to the sample, defined by the components of the 
highly anisotropic hole Lande tensor. We will assume that 
the system geometry is very close to the Voigt configura- 
tion, with the magnetic field only slightly tilted from the 
system plane, so that cos tp « and sin tp « 1. 

The system dynamics is a composition of many pro- 
cesses. For the optical hole-trion transition, we assume the 
radiative recombination rate 71 and the additional pure de- 
phasing rate 70. The hole and the trion spin precess in 
the magnetic field with the Larmor frequencies lo^ and 
cJt> respectively. The longitudinal relaxation rate for the 
hole spin is = n + + k_, where k_, k + are the rates 
for the spin flip transitions to the lower and to the upper 



Zeeman level, related by the detailed balance condition 
k„/k + = expfiujh/kBT), where k-g is the Boltzmann 
constant and T is the temperature. The hole spin transverse 
relaxation rate is = /2 + kq, where kq is the ad- 
ditional pure dephasing rate. We assume that the electron 
spin dephasing is slow compared to the trion lifetime so 
that the electron spin coherence is dominated by the latter. 

The simulations presented in this contribution are per- 
formed for T — 4 K and for the magnetic field of 7 T 
tilted at an angle ir/2 — t9 = 4° from the system plane. 
For the in-plane and perpendicular components of the hole 
g-tensor equal to 0.04 and 0.6, respectively, this yields the 
hole Larmor frequency us^ = 0.036 ps _1 and the hole spin 
quantization axis oriented at p = 44° from the normal di- 
rection. The trion Larmor frequency is 0.16 ps _1 (corre- 
sponding to the electron g-factor of 0.26). We will assume 
the decoherence times I/71 = 1.2 ns and T\ = 2 ns and 
no additional pure dephasing effects (70 = kq = 0). 

3 The FWM response. Our theoretical analysis is 
based on the method developed in Ref. [4j: the system evo- 
lution is studied in the density matrix formalism, with the 
optical and spin-related dephasing included via a Lindblad 
dissipator in the evolution equation. In the limit of coher- 
ent ultrafast optical driving, this approach yields analytical 
formulas for the single system dynamics, which allows one 
to perform averaging over an inhomogeneous distribution 
of various parameters in the ensemble. 

For a cr + probe, the only element of the density ma- 
trix just after the probe pulse that carries the e 2l< ^ 2 phase 
dependence is 



P3l(0+) 



p 13 (0-)e 2 *sin 2 ^, 



where (j>2 is the phase of the probe pulse, ot2 is its area, 
and 0~ denotes the time instant just before the arrival of 
the pulse. The evolution of the system state is then calcu- 
lated using the Master equation in the Lindblad form, like 
in Ref. Q. The er + interband coherence at a time t > is 

P3i(i) = ^rf±, ± e A ±^ t pi3(0 + )e~^, 
±,± 

where the amplitudes and the exponents are given by 
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Here E is the interband transition energy (at zero field) and 
the dephasing constants are 

r = 470 + 271 + k + k + + K- 
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and 



P = (k- — k+) [(cos tp — 1) sin 2 ip + l] . 



Now, we need to find pi3(0~). For a a + -polarized 
pump, just after the pump pulse (at t — — r + ), the only 
non-zero element linear in the pump amplitude is p^i. By 
solving the Master equation one finds at t = 0~ 



P13(0")=£ 



cue 



A*t iEr/h 



(1) 



where 

Ci = Pl3(-T + )di. 

Thus, at the time t one finds 

P31 (t) = isina 1 sin 2 ^e^-^) 



-iE(t-r)/h ' 



ij 



d t d ie xa+ ^ T . 



When averaging over the distribution of interband energies 
E in a typical QD ensemble, the exponent produces a nar- 
row echo peak around t = r. Since the spin-related evo- 
lution (frequencies A^) is slow on this time scale, we can 
put t = t under the summation. Then the magnitude of the 
time-integrated response is proportional to 

2 



FWM 



in 



COS U!tT)e 



-fx/2 



Ife 



-Pt/2 



(cos ip 

+ (cos<^ + l) 2 e /3r/2 + 2(1 - cos 2 ip)cosuj h T 



(2) 



In order to take into account also the inhomogeneity of 
the Lande tensors at different trapping centers this result 
should be averaged over a certain distribution of the Lar- 
mor frequencies oj^ and uj t , which we will assume to be 
Gaussian and characterized by the variances a\, of. We 
assume here that the distribution of Lande tensors is un- 
corrected to the spectral positions of the trion transitions. 

In Fig. 0I a X we show the calculated time-integrated 
FWM signal in the absence of inhomogeneous distribution 
of the Larmor frequencies. As follows from Eq. (0, oscil- 
lations at various frequencies are present in the optical sig- 
nal, corresponding to combinations of the hole and trion 
frequencies. A more transparent picture is obtained after 
performing a Fourier-transform of this signal, as shown in 
Fig-QIb). Here one can see the zero-frequency line, as well 
as lines at the frequencies ±wt, ±0Jh, and ±u> t ± ^h- 

In a real system, the optical response is affected by the 
inhomogeneous distribution of the relevant parameters. In 
the present case, the Larmor frequencies are usually not 
identical for each hole-trion system due to variations of 
the (^-factors in the nanostructure. In Fig. [2] we show the 
results of a simulation performed for the same parame- 
ters as in the previous case but with the additional effect 




Figure 1 The FWM signal from a system without inho- 
mogeneous broadening of the Larmor frequencies, (a) The 
time-integrated FWM signal as a function of the delay 
time; (b) The real part of the Fourier transform of the time- 
resolved signal. 
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Figure 2 The FWM signal from a system with inhomo- 
geneously broadened Larmor frequencies, (a) The time- 
integrated FWM signal as a function of the delay time; (b) 
The real part of the Fourier transform of the time-resolved 
signal. 



of inhomogeneous broadening of the Larmor frequency, 
which is assumed to be equal to 20% of their average val- 
ues. Due to this inhomogeneity effect, the oscillations in 
the time integrated response are damped and the signal is 
dominated by the monotonic decay (except for short delay 
times). Correspondingly, all the non-zero frequency peaks 
in the Fourier spectrum are strongly broadened. However, 
the central peak remains unaffected. 

4 Discussion and conclusions. An interesting 
point in the discussion presented above is the presence of a 
zero-frequency component which produces a central peak 
in the Fourier transform, composed of two Lorentzian con- 
tributions: one with the width (r + f3)/2 and another one 
with (r—f3)/2. Interestingly, as follows from Eq. ([TJ, there 
are no zero-frequency components in the polarization evo- 
lution between the pulses. Therefore, the non-oscillatory 
part of the FWM response can be interpreted as a result of 
partial refocusing of the Larmor precession by the probe 
pulse. 

The presence of this zero-frequency component, which 
is insensitive to the inhomogeneous distribution of the g- 
factor, opens a possibility of extracting useful information 
on the spin-related system parameters. If the trion life- 
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time is sufficiently long compared to the hole spin decoher- 
ence times (j3 not too small compared to F) then the two 
Lorentzian components of the central line, with the widths 
r ± j3 can be separated by fitting, from which the values 
of r and (5 can be deduced. The value of f3, along with 
the detailed balance relation, allows one to extract the rates 
K±, hence the longitudinal decoherence time. On the other 
hand, r involves both the spin-related rates K± t o and the 
optical lifetime and dephasing rates 70.1. The latter, how- 
ever, can often be deduced independently for a given sys- 
tem, which can allow one to find also the value of ko and 
to calculate the intrinsic transverse spin decoherence time 
T 2 . 

This should be compared with the information avail- 
able from a TRKR experiment |4|. There, the hole-related 
response originates from the occupations of the hole Zee- 
man sublevels after the trion recombination. Therefore, the 
TRKR experiment allows one to extract the hole spin- 
related information even if the exciton lifetime and co- 
herence time are very short, as was indeed the case in 
the experiments (HO ■ However, as the information on the 
transverse spin dephasing in that experiment is deduced 
from the spin precession signal, it is convoluted with the 
inhomogeneous effect which may even completely domi- 
nate the decay of coherent precession signal. Therefore, in 
a TRKR experiment essentially only the inhomogeneous 
transverse decoherence time T 2 * is available. 

Thus, we have shown that a FWM experiment per- 
formed on an ensemble of quantum dots or other trapping 
centers doped with excess holes can provide useful infor- 
mation on the properties of hole spin precession and de- 
coherence. Although this information is only available un- 
der favorable conditions related to the various decoherence 
rates in the system (long exciton coherence) it is insensi- 
tive to inhomogeneous effects, in particular to a variation 
of g-factors. Therefore, we conclude that the FWM exper- 
iments may be a useful tool to extract spin-related infor- 
mation which is complementary to that available from the 
TRKR study. 
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